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Abstract — Spatial sampling is traditionally studied in a static 
setting where static sensors scattered around space take mea- 
surements of the spatial field at their locations. In this paper 
we study the emerging paradigm of sampling and reconstructing 
spatial fields using sensors that move through space. We show 
that mobile sensing offers some unique advantages over static 
sensing in sensing time-invariant bandlimited spatial fields. Since 
a moving sensor encounters such a spatial field along its path as 
a time-domain signal, a time-domain anti-aliasing filter can be 
employed prior to sampling the signal received at the sensor. 
Such a filtering procedure, when used by a configuration of 
sensors moving at constant speeds along equispaced parallel 
lines, leads to a complete suppression of spatial aliasing in the 
direction of motion of the sensors. We analytically quantify the 
advantage of using such a sampling scheme over a static sampling 
scheme by computing the reduction in sampling noise due to 
the filter. We also analyze the effects of non-uniform sensor 
speeds on the reconstruction accuracy. Using simulation examples 
we demonstrate the advantages of mobile sampling over static 
sampling in practical problems. 

We extend our analysis to sampling and reconstruction 
schemes for monitoring time-varying bandlimited fields using 
mobile sensors. We demonstrate that in some situations we 
require a lower density of sensors when using a mobile sensing 
scheme instead of the conventional static sensing scheme. The 
exact advantage is quantified for a problem of sampling and 
reconstructing an audio field. 

I. Introduction 

The typical approach for measuring a spatial field makes 
use of static sensors distributed over the area of interest lU. 
Consequently much of the literature on spatial sampling and 
reconstruction have focused on such static sensing schemes 
II2II . I3I . An emerging paradigm in spatial sampling is the 
use of mobile sensors that move through the area of interest, 
taking measurements along their paths IH, O. Mobile sensing 
schemes have several advantages over static schemes, the chief 
of which is the fact that a single mobile sensor can be used to 
take measurements at several distinct positions in space. Such 
a scheme is often more cost-effective and easier to implement 
than static sensing since it requires only a single sensor for 
monitoring a large spatial area. In this paper we illustrate 
and analyze various unique aspects of mobile sensing and the 
advantages that they offer. 

Consider a time-invariant spatial fielcfll defined over d- 
dimensional space represented by a square integrable mapping 
/ : R with / G ^^(R^). For any r G R^, the quantity 

field that is a function of space alone and does not vary with time. 



/(r) represents the value of the field at spatial location r. 
The field could represent, for instance, some spatially varying 
parameter like the temperature of air or the concentration of a 
pollutant in the air. Suppose further that the field / is slowly 
varying in space and can hence be modeled as a spatially 
bandlimited field. Let / represent the observed field that is a 
noisy version of the field of interest expressed as 

/(r) = /(r)+^(r),r gR^ (1) 

where w denotes non-bandlimited spatial noise, which we 
refer to as environmental noise. The objective in a typical 
sampling and reconstructing scheme is to use the samples of 
the observed field / to obtain a reconstruction / of the field 
/ such that the mean-square error (MSE) E[|| / — / II2] in the 
reconstruction is minimal. If the noise w is absent, then the 
observed field / = / is bandlimited in space, and we know 
from classical sampling theory [6] that we can recover the 
field exactly from samples of the field taken by static sensors 
located on a lattice of points in space. However, iiw^{) then 
the observed field / is not bandlimited and we expect to see 
some effects of spatial aliasing while sampling the field using 
static sensors. More importantly, unlike in the case of sampling 
a time-domain signal, there is no way to implement a spatial 
anti-aliasing filter in a static sensing setup. This drawback 
of spatial sensing using static sensors has been observed in 
various applications (see, e.g., Q, ID, 191). 

Sampling using mobile sensors provides us with an ap- 
proach that partially addresses the issue of spatial aliasing. 
Mobile sensing allows filtering in time prior to sampling which 
induces filtering over space in the direction of motion of the 
sensor. Such spatial filtering is not possible in a static sensing 
setup. A moving sensor receives as input a time-domain signal 
representing the field along the path of the sensor given by, 

m = f{r{t))=f{r{t))+w{r{t)) (2) 

where r{t) G R^ denotes the position of the sensor at time t. 
The analog signal s{t) can be passed through an analog anti- 
aliasing filter prior to discretizing into samples. Such filtering 
discards out-of-band noise in the direction of motion of the 
sensor thus inducing spatial smoothing. Implementing such a 
filter requires redesigning of the sensing process employed in 
the sensor and may be feasible only in some scenarios. For 
the purpose of illustration consider a problem of measuring 
the concentration of a gas in the air along a straight road, 
assumed to be constant in time. Assume further that one has a 
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sensor that can measure the average concentration of the gas 
in a chamber. We want to design a scheme that computes the 
filtered samples 

/CO 
s{T)h{nT — T)dT 
-oo 

at times nT for all n G Z. If we assume that is a finite 
impulse response filter satisfying h{t) > for all t and 
h{t) ={) for \t\ > ^ , then such a filter can be implemented 
by designing a system in which a chamber is mounted on 
a vehicle moving along the road at a constant velocity. Air 
is pumped into the chamber through an opening whose size 
can be varied dynamically such that the rate of entry of air 
at time r is proportional to h{T^^ + — r) where \_x\ 
indicates integer component of a real number x. The chamber 
is periodically emptied after each sample is measured at times 
(n+ |)T for all n e Z. Such a scheme essentially implements 
an analog domain sampling prefilter with impulse response h 
up to a normalization constant. This scheme can be generalized 
to more general non-negative finite impulse responses h with 
wider support, if we use multiple chambers. 

However, in the proposed approach, a caveat to note is 
the peculiar fact that such filtering permits spatial smoothing 
only in the direction of motion of the sensor. Hence, for 
spatial fields of dimension d > 1, the anti-aliasing filters are 
thin, i.e., the effective spatial impulse response of the filter is 
supported on a set of dimension 1. Thus this form of spatial 
smoothing allows us to discard only one component of the 
out-of-band noise. For the sake of illustrating the potential 
advantage offered by such a filtering scheme, let us consider 
the problem of sampling a field in two-dimensional space 
having a Fourier transform that is bandlimited only in one 
direction as shown in Figure |l(a)[ Sampling such a field 
using static sensors will always lead to aliasing because the 
repetitions in the sampled spectra necessarily overlap. An 
example of the sampled spectrum obtained by static sampling 
on a lattice of the form {{mAx^nAy) : m,n G Z} is shown 
in Figure [T(b)| where Ay < However, we will see later that 
such a field can be sampled on the same lattice using sensors 
that move along equispaced straight lines parallel to the x-axis 
as shown in Figure [2(a)| and use ideal anti-aliasing filters in the 
time domain. This leads to a complete suppression of aliasing 
as shown in Figure |l(c)| We analyze such a mobile sensing 
scheme for sampling bandlimited fields later in the paper and 
quantify advantages obtained in terms of suppressing out-of- 
band noise. 

The scenario is a little different in the case of sampling 
time- varying bandlimited fields, represented by functions of 
both time and space of the form /(r, t) where r denotes 
position and t time. Here the advantages of mobile sensing 
are less pronounced since the field values at various points 
in space are also varying in time. Hence it is possible to 
filter in time even with static sensors. Furthermore, one has 
to account for the Doppler effect while sensing with moving 
sensors (see e.g., |[TOl sec 5.2]). Nevertheless we show that in 
some scenarios mobile sensing requires fewer sensors than the 
number of sensors required with static sampling. These results 
could have potential applications in limiting spatial aliasing in 
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TABLE I 

Summary of model assumptions in different sections. The last 
two columns represent two different additive noise models 
that we consider - analog spatial noise added to the field 
prior to sampling, and discrete noise added to the field 
measurements obtained after sampling. 



certain sampling problems, e.g., wave-field synthesis lO. 

The advantage of mobile sensing over static sensing in 
suppressing spatial aliasing has been noted in the context of 
audio source localization in where the authors show that 
using a planar rotating array of microphones, the effective 
number of measurement points can be increased thus reducing 
the amount of spatial aliasing while using the beamforming 
technique. In some other works moving microphones have 
been used for estimating room impulse responses 111] and 
head related impulse responses fT2l . In these works the desired 
responses are estimated from a finely sampled version of the 
audio signal received at the moving microphones. Other works 
on mobile sensing focus on adaptive path-planning algorithms 
El, ifTSl for environmental monitoring. However, to the best of 
our knowledge, this is the first work to illustrate the possibility 
of implementing spatial smoothing by using a mobile sensing 
scheme together with time-domain anti-aliasing filtering, and 
to quantify the improvements of such a scheme over static 
sensing. In earlier work 1141, we studied the problem of 
designing trajectories for mobile sensing which minimize the 
total distance required to be traveled by the sensors per unit 
area of the field being sampled. 

The paper is organized as follows. We study the sampling 
of time-invariant fields in Section [III and time- varying fields 
in Section [Till We discuss sensor trajectories, filter designs, 
reconstruction schemes, and comparisons with static sensing 
in various aspects. In order to highlight various advantages of 
mobile sensing we have considered several different models 
for the field and the noise. To enable easy navigation through 
the paper we summarize the model assumptions used in the 
various subsections in Table U In Section [IVl we discuss a 
simulation example comparing mobile and static sampling for 
a practical problem and conclude in Section |Vl 

II. Time-invariant fields 

Consider a time-invariant field in d-dimensional space rep- 
resented by a mapping / : ^ R with / G ^^(R^). We 
define its Fourier transform F as 

F{uj)= /(r)exp(-i(a;,r))dr, a; G R^ 

where i denotes the imaginary unit, and (u^v) denotes the 
Euclidean inner product between vectors u and v. We assume 
that the field / is bandlimited to a set ^7 C R^, i.e., suppose 
that the Fourier transform F of / is supported on a known 
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(a) Spectrum of field bandlimited to M x [—p,p]. (b) Aliased spectrum under static sampling. (c) No aliasing under mobile sampling with fil- 
tering in the x-direction. 



Fig. 1. Sampling a two-dimensional field bandlimited only in one direction. Spectrum under static sampling is aliased but that under mobile sampling is not. 



set n C M^, so that F{uj) = for uj ^ Q. We use / to 
denote the observed field, which may either be the field / 
itself or a noisy version of it as shown in ([T]). In the noisy 
case, we model the environmental noise w{r) as a zero mean 
wide sense stationary (WSS) process with unknown power 
spectral density Sw For mathematical regularity, we assume 
that the noise power spectral density decays at a faster rate 
than 0(||cj||^^), i.e., we assume that 

5^(0;) = for some r > d. (3) 

Time-invariant fields are particularly well- suited for mobile 
sensing schemes since the field does not vary as the sensor 
moves around taking measurements in space. We distinguish 
between two distinct scenarios - one-dimensional fields in 
which a single moving sensor can visit all points in the one- 
dimensional spatial region of interest, and higher-dimensional 
fields in which each sensor can measure the field only on a 
one-dimensional sub-manifold of the spatial region of interest. 
To illustrate this difference we consider one-dimensional fields 
and two-dimensional fields with the understanding that the 
analysis for two-dimensional fields can be easily generalized 
to higher-dimensional fields. In this section, we detail the 
mobile sensing scheme and describe how a time domain anti- 
aliasing filter can be used to perform spatial smoothing. We 
quantify the improvements of such a sensing scheme over the 
static sensing scheme. We initially assume that the sensors are 
moving with constant velocities and later consider the scenario 
with non-uniform speeds. 

A. Sensor trajectories for mobile sensing 

For sampling a field in where d > 2, there are several 
possible choices of trajectories that can be used by the moving 
sensors. In our recent work UlAj we studied the problem of 
designing sensor trajectories that admit perfect reconstruction 
of bandlimited fields from measurements taken by the sensors 
moving along these trajectories in a noise-free setting. We 
introduced the notion of optimal trajectories that minimize the 
total distance required to be traveled by the moving sensors. 
For d = 2 we showed that a set of trajectories comprising one 
set of equispaced parallel lines is optimal from certain classes 
of trajectories. In this paper, we will restrict ourselves to such 
a collection of trajectories for fields in and assume that we 



have one sensor moving along each line taking measurements 
on its path. 

B. Sampling and reconstruction 

Consider a sensor moving at a constant velocity through 
space. The position r{t) of the sensor at time t is given by 
an affine function of the form r{t) = u -\- vt where u^v G 
R^ represent the initial position and velocity of the sensor 
respectively. Without loss of generality we assume that u = 
for simplifying the analysis. In the absence of noise the time- 
domain signal seen by such a sensor is given by 

so{t) = f{r{t))=f{u^vt),teR. 

In this scenario it easily follows (see, e.g., Lemma 2.2]) 
that the signal so(.) is bandlimited to 

Qso '= {{v^oj) :uj en} CR. (4) 

It follows via the Nyquist sampling theorem that sq{.) can 
be perfectly recovered by sampling it uniformly at tempo- 
ral intervals less than or equal to 27r/(maxl]sQ — minl^gQ). 
Furthermore, in the presence of noise, the signal received by 
the sensor can be passed through an anti-aliasing filter with 
passband aligned with prior to sampling. This limits the 
contribution of out-of-band noise in the samples. 

1) One -dimensional field: Suppose /(.) denotes a one- 
dimensional field bandlimited to 1] = [— p, p]. In this case, if a 
sensor moves along the field at a constant velocity v, the signal 
it sees in the absence of noise is bandlimited to [—vp^vp] as 
we argued in (|4]). Hence, in the presence of noise, the signal 
can be filtered prior to sampling. Let h{.) denote the impulse 
response of an ideal filter with passband in the interval [— p, p] : 

h{x) = -sinc( — ), X eR (5) 

TT TT 

where sinc(x) := ^^"^^ . Then the ideal choice for an anti- 
aliasing filter while sampling the time-domain signal s{t) = 
f{vt) is given by /laa(^) = vh{vt). Let s{t) denote the signal 
at the output of the anti-aliasing filter. Suppose the sensor takes 
measurements every T time units after passing the observed 
signal through the filter. Thus we get uniform samples of 

sit) = (s * /laa)(^) = (^0 * h^^)(t) iw ^ h^^)it) 
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where w{t) = w{vt) and * denotes convolution. Since / is 
bandlimited we can write 

s{t) = f{vt)^{w^h){vt). 

Now, if the sensor takes samples every T time units, we 
essentially get uniform noisy samples of the field at intervals 
of vT spatial units. We know from classical sampling theory 
ifTSll that when noise is absent the field /(.) can be exactly 
recovered from these samples provided that the sampling 
interval satisfies T < Furthermore, in the noisy case, out- 
of-band noise can be suppressed in the reconstruction of the 
field by using sine interpolatiorH 



fix) = y: —4jt) 

jez 



smc 



p{x - jvT) 



X eR (6) 



provided T < This interpolation is well-defined when the 
noise satisfies ([3]). However sine interpolation is sensitive to 
errors in the samples and can lead to unbounded errors in the 
approximation at some values of x. This can be avoided by 
using alternate kernels in place of the sine kernel (see, e.g., 
ifTvl and references therein). 

2 ) Two-dimensional field: In the case of a two-dimensional 
field /(.), we consider sensors moving along equispaced par- 
allel lines through space. For a given Vt the optimal orientation 
of these parallel lines can be computed as described in llT4l . 
However, to simplify analysis we assume that the sensors 
are moving at a constant velocity v along lines parallel to 
the X-axis spaced Ay units apart as illustrated in Figure 
|2(a)| The position of the j-th sensor at time t is given by 
{vt^jAy). Hence the j-th sensor is exposed to the signal 
Sj (t) = f{vt^jAy). In the absence of noise, we know from dU 
that these signals are bandlimited. Hence in the noisy scenario 
the signals sj {t) can be filtered and sampled uniformly, just 
like in the case of the one-dimensional field as we outlined in 
Section ITl-Bll In particular, if ft = [— p, p] x [— p, p], then 
it follows from ^ that the signal so{t) = f{vt^jAy) is 
bandlimited to an interval of the form [—pv^ pv]. Hence, like in 
Section Hl-Bll the filter /laa(^) = vh{vt) can be used prior to 
uniform sampling. The frequency response of this ideal time- 
domain filter is shown in Figure |2(b)| Thus, at the output of 
the sampler we obtain uniform samples of 

Sj{t) = f{vtJAy) ^ {wj ^ h){vt) 

where Wj{x) = w{x^jAy) and h is the filter defined in ©. 
Suppose that the samples are taken at time-intervals of T < 
— units. Let A^^ = vT. We further assume that the samples 
taken by the sensors are aligned with each other, and hence 
the collection of samples from all the sensors lie on a two- 
dimensional lattice of the form {{iAx^jAy) : z, j G Z} and 
can be expressed as 

Sj{iT) = f{iAxJAy) + {wj * h){iAx), ij e Z. (7) 

^Note that our rationale behind using sine filters for anti-aliasing and 
interpolation is the fact that the only assumption we have about the field is 
that it is bandlimited. For a stochastic field model the MSB can be minimized 
by using a Wiener filter fT6l . 
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Fig. 3. Sampling and reconstruction setup. 



Thus, results from classical sampling theory f6l can be used 
to estimate the field from these samples as: 



AxAyp'^Sj{iT) . (p{x-iAx) 
^ — ^smc ' ^ 



smc,'^^(^^i^~) (8) 



provided Ax, Ay < ^. The reconstruction given in ([S]) is well- 
defined when the noise satisfies © and is exact when noise 
is absent. 

Using the notation /i[n] = Sj{iT) where n = the 
filtering and sampling relation of d?]) can be written as a two- 
dimensional convolution as follows: 

l^[n\ = (/^/is)(An),n G Z^, 

where ^ denotes two-dimensional convolution, 

A., 



A 



A,. 



and hs represents the effective two-dimensional sampling 
kernel induced by the sampling trajectories of Figure [2(a)| and 
the filtering operation of Figure |2(b)| given by 



hs{x,y) = -sine ( — ) S{y) 



(9) 



where 5{.) represents the Dirac delta function. This kernel 
has the following representation in the Fourier domain as 
illustrated in Figure |2(c)| 



_ j 1 for < \ujx\ < p 
) else. 



(10) 



Figure [3] shows the system of sampling and reconstruction 
where u plays the role of / and z> the role of /. In the figure the 
samples /i[n] passes through a discrete to continuous converter 
to produce /is(^), an impulse stream in the continuous space 
given by 

/is(r) = /i[n]J^(r — An) 

where 6'^ represents the Dirac delta function in two dimen- 
sions. The reconstruction kernel hr{.) is the two-dimensional 
sinc-kernel 

Mx..) = ^sinc(5i)si„c(^). 

The kernel has the following representation in the Fourier 
domain 



AxAy for < \ujx\, \oJy\ < p 



r ^x^y ... 

^ ^ 1^ else 



(11) 
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Hs{u) = 



V (velocity) 

—vp vp 

(a) Equispaced parallel line trajectories, (b) Frequency response of time domain anti- (c) Frequency response of the induced sam- 

aliasing filter. pling kernel is supported on {uu : \ujx\ < p}- 



Fig. 2. Sampling a two-dimensional field using mobile sensors: Sensor trajectories, frequency response of time-domain filter, and the induced sampling 
kernel. 



where u = (ux^oOy). In this representation the reconstructed 
field of ^ is given by 

/(r) = Yl /^W^r(r - An),r G 

We will use this new representation for simplifying the dis- 
cussion in the rest of the paper. 

C. Comparison with static sampling: Aliasing suppression for 
non-bandlimited fields 

Before proceeding to discuss bandlimited fields in detail, 
we now take a slight detour to consider the problem of 
sampling a non-bandlimited field in a noise-free setting. We 
know from classical sampling theory that sampling a signal 
at a rate less than the Nyquist rate leads to aliasing in the 
reconstructed signal, which is highly undesirable. In practice, 
while sampling a non-bandlimited signal in the time-domain, 
one typically employs an anti-aliasing filter to suppress the 
out-of-band portion of the signal. The ideal choice for the anti- 
aliasing filter and the reconstruction filter are ideal low-pass 
filters with cutoff frequencies given by half of the sampling 
rate. However, for sampling in the spatial domain, it is not 
possible to implement a spatial anti-aliasing filter in a static 
sensing setup. Nevertheless, in a mobile sampling scheme, an 
ideal time-domain filter can be used to reduce the amount of 
aliasing. In the one-dimensional case, as we argued in Section 
III-Bll this sampling procedure effectively implements an ideal 
anti-aliasing filter thereby suppressing aliasing completely. 
Thus the mobile sampling scheme is completely devoid of 
aliasing effects even while sampling a non-bandlimited field. 
Furthermore, the squared error in the reconstruction can be 
further reduced by simultaneously increasing the sampling rate 
and the bandwidth of the anti-aliasing filter. 

The scenario is different in the two-dimensional case. Sup- 
pose we design the sampling scheme for fields bandlimited 
to ^7 = [—pjp] X [— p, p] and suppose that the observed 
noise-free field / = / is in fact bandlimited to the bigger 
sell ^ = [—p^p] X [—p^p] where p > p. Such a field is 
shown in Figure |4(a)[ In a static sensing scheme this field 

^The assumption that the field is bandlimited to Cl is only used to simplify 
the illustration. The advantages of filtering persist even when / is not 
bandlimited. 




□ Spectrum of interest 

□ Out-of-band energy 

(a) Non-bandlimited spectrum. 



□ Spectrum of interest 

□ Out-of-band energy 

(b) Filtered spectrum. 



Fig. 4. Spectrum of a non-bandlimited two-dimensional field and its filtered 
version. 



is sampled without filtering and hence the sampled spectrum 
and reconstructed spectrum are aliased in both the ujx and ujy 
directions as shown in Figures |5(a)| and |5(b)| In the mobile 
sampling case, however, we can use the anti-aliasing filter 
described in Section III-B2I In this case we know from the 
form of the effective sampling kernel in Figure |2(c)| that 
the out-of-band energy in the ujx direction is filtered out 
as shown in Figure |4(b)| Hence the sampled spectrum and 
reconstructed spectrum are aliased only in the ujy direction 
as shown in Figures |5(c)| and |5(d)| respectively. Although the 
obtained reconstruction / is not completely devoid of aliasing, 
the squared error ||/ — /II2 in the reconstruction is significantly 
lower than the squared error in the reconstruction obtained 
via static sampling which is aliased both in the ujx and ujy 
directions. 

Two observations are in order. Firstly, from the discussion 
above and Figure [5] it is clear that the scheme of mobile 
sampling with filtering is effective in suppressing aliasing 
in the cUx direction even when the field is not bandlimited 
in the uJx direction. We discussed such an example in the 
introduction. The Fourier transform was supported on a semi- 
infinite set and the sampled spectrum is completely devoid of 
aliasing as illustrated in Figure [T] Such a complete suppression 
of aliasing is not possible with static sensing unless the field 
is bandlimited to a compact set in . Secondly, it is possible 
to improve upon the mobile sampling scheme described here 
by employing mobile sampling along more complex sampling 
trajectories. For instance, in the example of Figure H] if we 
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replaced the sensors moving parallel to the x-axis with a 
similar set of sensors moving at constant velocities parallel 
to the y-Sixis then it would be possible to suppress aliasing in 
the ujy direction by using ideal anti-aliasing filters as before. 
Now if we had both kinds of sensors, those moving along the 
X direction and those moving along the y direction, then the 
prefiltered samples from the former set of sensors are devoid 
of aliasing in the ujx direction and the prefiltered samples 
from the latter set of sensors are devoid of aliasing in the 
ujy direction. The samples from both these sets of sensors can 
be combined to obtain better aliasing suppression than that 
obtained in the Figure |5(d)| A simple way to do this would 
be to first obtain two different reconstructions of the field, the 
first using the set of samples from the sensors moving along 
the X direction and the other from the sensors moving along 
the y direction. Now the aliased portions of each of these 
reconstructions can be filtered out and the resultant fields can 
be averaged to obtain a reconstruction with less aliasing than 
either original reconstruction. 

D. Comparison with static sampling: Noise suppression for 
bandlimited fields in noise 

We now study the reduction of out-of-band noise that can 
be obtained in the mobile sensing scheme using an anti- 
aliasing filter while sampling a bandlimited field in WSS 
noise satisfying (O. The analysis is straightforward for one- 
dimensional fields. We provide a more detailed study for the 
two-dimensional case. 

1 ) One-dimensional field: In a scheme of sampling a one- 
dimensional field using static sensors, one has access only to 
field measurements taken on a discrete set of points. Hence it 
is not possible to completely filter out the out-of-band noise in 
a static setup. However, if we use a sensor moving at constant 
velocity it is possible to filter out all the out-of-band noise as 
we showed in Section lll-Bll We can quantify the advantage by 
comparing the noise variance in the reconstruction of (O with 
the noise variance in the reconstruction of a static sampling 
scheme. Assume that field / is bandlimited to 1^ = [— p, p] and 
is sampled at the Nyquist rate. In the case of static sampling 
there is no anti-aliasing filter and hence the noise variance in 
the resulting reconstruction is given by 

1 

where Su,{') denotes the power spectral density (p.s.d.) of the 
noise process w{.). In the mobile sampling case, however, the 
field is effectively filtered by an anti-aliasing filter given by © 
as in Section III-Bll prior to sampling. Since the anti-aliasing 
filter prevents aliasing due to out-of-band noise, the total noise 
in the reconstruction of (O is given by 



1 



(13) 



Clearly, we see that the contribution from the repetitions of S^^ 
that appear in the static reconstruction of ([T2l) is absent in the 
variance of noise in the reconstruction of ([T3]) corresponding 
to the mobile sampling and filtering scheme of Section III-Bll 



Moreover, the contribution from the terms with n 7^ in ([T2l) 
leads to aliasing in the reconstruction which is particularly 
undesirable. 

2 ) Two-dimensional field: Consider a two-dimensional spa- 
tial field / bandlimited to Q = [-p,p] x [-p,p] C M^. 
The static and mobile schemes for sampling the field can 
be represented as shown in Figure [3] where u represents the 
input field being sampled, hs denotes the specific choice of 
the sampling kernel and A is a 2 x 2 matrix that generates the 
spatial sampling lattice. In general, u is 3. noisy version / of 
the field /. 

We consider a static sampling scheme that uses a rectangular 
sampling grid at the Nyquist sampling rate i.e., = = ^ 
and we consider a mobile sampling scheme as in Section llI-B2l 
with the spacing A^ between the sensor trajectories equal to 
the Nyquist interval ^. We consider two extreme choices for 
the low-pass filter employed prior to sampling in the mobile 
sensing scheme - the ideal low-pass filter (LPF) with a sine 
response, and a more practically feasible filter, the box-filter 
whose impulse response is a rectangular pulse. 

We first characterize the response of the sampling schemes 
to noise. Suppose the input u in Figure [3] represents weak- 
sense stationary noise with power spectral density Sj,(cj),cj G 
R^. Then the p.s.d. of p is given by. 



nG2 



S^(A-i(cj - 27rn))|i^s(A"H^ - 27rn))|2 
A^A^ 



The output signal z> is cyclo stationary since its autocorrelation 
function is invariant to shifts by 27rA~^. Its effective p.s.d. 
can be computed just as one computes the power spectra of 
linearly modulated signals (see, e.g., (TSl Sec. 3.4-2]) and is 
given by 



A^A, 



(14) 



Below, we separately compute the noise variances in the case 
of static and the two cases of mobile sampling. 

Static sampling 

In the case of static sensing we do not have any spatial 
filtering and hence the sampling kernel is an impulse function 
hs{x^y) = S{x)S{y) and is equivalently given in the Fourier 
domain by, 

Hs{uj) = 1 for all a; G (15) 

Since we assume Nyquist rate sampling we also have A^^ = 

Ay = -. Substituting ([T5]) in ([T4l) and integrating we obtain 

^ ^2 



the following expression for the noise variance cr^^^^. 
(^Itat = J-2 [ XI S^(cj - 2np)duj. 



(16) 



Since we are sampling at the Nyquist rate we do not observe 
any signal distortion in this case. 

Mobile sampling with ideal low-pass filter 

If the ideal sine filter is employed as the LPF prior to sampling 
at the mobile sensors, the effective two-dimensional sampling 



□ Non-overlapped spectrum 
1} Overlapped portions 



□ Unaliased portion 
Ij^ Aliased portions 



□ Non-overlapped spectrum 
■ Overlapped portions 



□ Unaliased portion 
■ Aliased portions 



(a) Sampled spectrum from static (b) Reconstructed spectrum from (c) Sampled spectrum from mobile (d) Reconstructed spectrum from 
sampling. static sampling. sampling with filtering. mobile sampling with filtering. 



Fig. 5. Comparison of aliasing in static and mobile sampling schemes while sampling a two-dimensional field. The reconstructed field is aliased in two 
directions for the static scheme while only in one direction for the mobile scheme. 



kernel hs is given by (O which is equivalently given in the 
Fourier domain by, 



Hs{uj) 



1 forO<|a;^|<p 
else. 



(17) 



We note that the filter Hs is frequency limited in the x- 
direction. Further, since the spacing between the trajectories is 
equal to the Nyquist interval A^y = ^ it follows that as long as 
the sampling interval along the trajectories satisfies Ax < ^, 
we have the following expression for the noise variance a^^: 



"^^^ nez 



[UJ — 



(18) 



Comparing the expressions ([T6l) and (fTSl) we see that in the 
case of mobile sampling, only spectral shifts at lattice points 
along the x-axis contribute to the noise variance whereas in 
the case of static sampling spectral shifts from all points 
in the two-dimensional lattice contribute to the noise signal 
spectrum. The exact value of the reduction in noise variance 
obtained with mobile sampling can be computed if the true 
noise spectrum Su{-) is known. 

Mobile sampling with box filter 

We now consider the case of mobile sampling with a box 
filter that is easier to implement. In this case we expect some 
distortion in the reconstruction because the filter response 
is not flat in the pass-band. The effective two-dimensional 
sampling kernel in this case is hs{x^ y) = bo^/\^^{x)5{y) where 



boxA,(x) := 







\x\ < 
else 



where 2 A5 denotes the spatial width of the box-filter response 

and hi = (2p)i X!^ sinc^ {^) duj] \ The sampling 
kernel is equivalently given in the Fourier domain by 

K.Ab . AbUJx 
Hs{uJ) = smc . (19) 

TT TT 

The choice of ensures that Hs{uj)'^ duj = which is 
consistent with the responses in ([T5]) and ([TT]) . 

Since the box filter is a non-ideal LPF, we expect some 
distortion in the reconstructed field. As before, let us assume 



that Ax < - and that = ^. Let F denote the Fourier 
transform of the bandlimited field being sampled and F denote 
the Fourier transform of the reconstruction obtained by using 
an ideal reconstruction filter of the form (fTTI) . We have 

F{co) = ^F{lo) sinc(^^), co & Q. (20) 

TT TT 

Clearly, we see that the reconstructed field is a distorted 
version of the original field even in the absence of noise. The 
amount of distortion introduced for a given F can be computed 
using relation (|2Q|) . We can also quantify the variance cr^2 
the noise in the reconstruction using relations ([T9l) and ([T4l) : 



^m2 



47r4 



smc ( — [Ux — 

TT 



2iTnx 



))duj. 



Now if it is inexpensive to increase the sampling rate used by 
each moving sensor, we can let A^^ 0, whence we obtain. 



2 

^m2 



47r4 



/ y S,(w - (0, 2np)^)sinc2(^ 



)duo. 
(21) 



If we now also allow A 5 0, it is easy to see that 
1. Hence the reconstruction in (|2Ql) becomes accurate 
and the expression in ([2T1) reduces to that in ([TSl) . This 
means that we have limAb^o F{uj) = F{uj) for G ft, 
and limAb^o limA^^o crm2 = ^mi- This means that if the 
sensors oversample at high rates along their path and use a 
rectangular LPF with a short impulse response, we can recreate 
the performance with an ideal LPF and Nyquist sampling. 

We note that if the noise spectral density were white (i.e. 
Szy(cj) = 1 for all UJ G M?), then the expressions for the 
variance in ([T6l) , (fTSl) and (|2T1) tend to infinity. This is because 
unfiltered white noise samples have infinite variance. However, 
in practice environmental noise is never completely white. 
If the noise spectrum is flat with a bandwidth along each 
dimension equal to a times the field bandwidth, then we obtain 
a reduction in noise variance by a factor of a when we employ 
mobile sampling in place of static sampling, as shown below. 
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Proposition 2.1: Suppose the noise spectral density S^i-) 
takes a value of unity on the set [— ap, ap] x [— ap, ap] and 
elsewhere for some p > 0. Then the variances in ([T6b and ([TSl) 
decay as O(a^) and 0(a) respectively. In particular when a 
is an odd number the variances reduce to 



p^a^ 



and 



a: 



ml 



p'^a 



Thus if the noise spectrum is flat with a bandwidth along each 
dimension equal to a times the field bandwidth, then we obtain 
a reduction in noise variance by a factor of a when we employ 
mobile sampling in place of static sampling. 

E. Noise suppression via oversampling on sensor paths 

A unique aspect of a mobile sensing scheme is the fact 
that it is possible to sample at high rates along the paths of 
the mobile sensors. This is a significant advantage over static 
sensing, because in static sensing the sampling rate can be 
increased only by increasing the number of sensors which 
is, in general, more expensive than increase the density of 
samples taken by a moving sensor. Such oversampling helps 
in suppressing measurement noise, i.e., noise added to the 
discrete measurements after sampling. In the case of sampling 
a one-dimensional field bandlimited to (7 = [— p, p], as we 
discussed earlier, the Nyquist sampling rate is ^ where v is 
the velocity of the sensor. It is known (191 P- 136] that, in 
the presence of additive zero-mean white measurement noise, 
oversampling by a factor of k relative to the Nyquist rate 
leads to a reduction of noise variance in the reconstruction 
by a factor of 1/k. Similarly, with non-linear processing 
quantization noise can be reduced by a factor of (see, 
e.g., II20II , II2TII . and references therein). 

A similar noise reduction by oversampling can also be ob- 
tained while sampling a two-dimensional field. As in Section 
III-D2I let = [— p, p] X [— p, p] and assume that the sensors 
move along equispaced straight lines parallel to the x-axis 
and spaced A^y = ^ apart. The temporal Nyquist sampling 
rate along the lines is ^ where v denotes the speed of the 
sensors. Now if the sensors take samples at k times the Nyquist 
rate, the effective sampling operation can be interpreted as a 
measurement of the inner products of the field with vectors 
from k disjoint orthogonal bases, as discussed in |[T9l p. 136]. 
Since these vectors form a tight frame with redundancy k, it 
follows via |[T9j Prop. 5.3] that oversampling by a factor of 
k leads to a reduction in noise variance in the reconstruction 
by a factor of 1/k when sampling in the presence of additive 
zero-mean white measurement noise. It is possible that the 
known result 1221 on the reduction of quantization noise by 
oversampling can also be extended to the two-dimensional 
case, but we do not consider such an extension here. 

F. Non-uniform sensor speeds 

In practice it may not always be possible to ensure that 
the sensors move at a constant velocity. We now analyze the 
mobile sensing of a one-dimensional field / bandlimited to 
Vt = [— p, p] using a sensor with time-varying speed. The 
analysis can be generalized to sensors moving along straight 



line trajectories in higher-dimensional fields, since the restric- 
tion of the higher-dimensional field to such trajectories is also 
bandlimited (see, e.g., 1231 ). Suppose that the position x{t) 
of the sensor at time t is a known monotonically increasing 
function of time t. We have so far studied the case when x{t) 
is an affine function, i.e., when the sensor has constant speed. 
In that case the signal 5o(t) = f{x{t)) is exactly bandlimited 
and hence can be sampled uniformly. For non-affine functions 
x{t), the signal SQ{t) = f{x{t)) is not bandlimited. However, 
such a signal is also essentially bandlimited to some po < 00. 
The value of po can be calculated as we show in the appendix. 
Hence we can obtain an approximate reconstruction of so{t) 
by uniformly sampling an appropriately filtered version of it. 
Although this filtering operation is linear and time-invariant, 
there is still some distortion introduced because of the non- 
linearity in x{t) which we now analyze. Let hip{.) denote the 
impulse response of an ideal low-pass filter (LPF) employed 
prior to sampling. Let Aip denote its 3-dB spread in the 
temporal domain. Denote by z{t) the output of the low-pass 
filter. We have 



m 



s{T)hip{t — r)dr = z{t) + w{t) 



(22) 



with 



z{t) 



-L 



j f{x{T))hlp{t-T)dT 

" fix)hipit - T{x)) 



dx 



(23) 



v{T{x)) 

where v{t) = is the velocity function, T(.) := x~^{.) 

is the inverse function of x(.), and w{t) = w{x{r))hip{t — 
r)dr. Thus at the output of the sampler, we get uniform 
samples {zn := z{tn)} at times tn := nT. Two observations 
are in order. The spatial separation between two successive 
samples is approximately, x^+i — Xn ^ v{tn){tn-\-i — tn) 
which means that samples are taken farther apart in space 
when the sensor is moving fast. We also note from (l23l) that the 
effective spatial spread of the sampling kernel while obtaining 
sample Zn is given by v{tn)Aip which means that samples 
obtained while the sensor is moving fast are obtained via a 
broader effective sampling kernel in space. It is also clear from 
([23l) that the sampling kernel is also scaled down by a factor 
proportional to the velocity. Now, since z{t) is bandlimited, 
it can be reconstructed exactly from its uniformly spaced 
samples by sine interpolation. Hence we can reconstruct an 
estimate for the field as f{x) := z{T{x)). Clearly f{x) can 
be expressed as f{x) = z{T{x)) + w{T{x)) where the first 
term z{T{x)) represents the contribution of the true field in 
the estimate and the second term represents the contribution 
of noise. We note that even in the absence of noise the 
reconstructed field is a distorted version of the field due to 
the non-linearity in x{t). The distortion in z{T{x)) can be 
quantified as follows: 

/ (fix) - z{Tix))fdx = f (soit) - z{t))Mt)dt 

< v\\so-z\\l (24) 

where v = sup^ \ v{t)\ denotes the maximum speed of the 
sensor and ||so — ^||2 denotes the total energy in so{t) outside 
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of the passband of the LPF. This suggests that the amount of 
distortion in the reconstruction can be reduced by increasing 
the bandwidth of the low-pass filter. This however comes at 
a cost of increasing the contribution of noise w{.) in the 
reconstructed field /(.). As a heuristic one can use the effective 
bandwidth of the signal so{t) = f{x{t)) as the bandwidth of 
the low-pass filter. 

III. Time-varying fields 

We now consider the more general problem of sampling 
time- varying spatial fields using mobile sensors. We focus on 
time- varying fields in one-dimensional space - i.e., fields of 
the form f{x,t) where x G M is a one-dimensional spatial 
parameter and t denotes time. Our approach can be extended 
to time- varying fields in higher-dimensional spaces. 

A. Sampling and reconstruction 

Let f{x^t) where G M, denote a time-varying field in 
one-dimensional space. Suppose that / is bandlimited to 1] C 
M?. If x{t) denotes the position of a moving sensor at time 
t, the signal seen by the sensor is given by f{x(t)^t). For a 
sensor moving at a constant speed v the position is an affine 
function of time of the form x(t) := u ^ vt. Then it is clear 
from dll) that the signal so(t) = f{x{t)^t) is bandlimited to 



(25) 



Thus in the absence of noise, the signal so{t) is bandlimited 
and can be exactly reconstructed by taking its samples at 
uniform intervals, as in the time-invariant case we considered 
in Section III-B2I In the presence of noise, an anti-aliasing 
filter with the appropriate bandwidth can be employed prior 
to sampling like in the time-invariant case. Furthermore, if the 
sensors are moving at non-uniform speeds then the signals 
are not exactly bandlimited but they can be approximated by 
bandlimited signals by following an approach similar to that 
in Section Hm 

We consider a scheme of sampling using a uniform col- 
lection of mobile sensors moving with equal velocities and 
separated by a constant separation in space. Such a configura- 
tion of moving sensors is illustrated in Figure |6(a)| Each line 
in the figure represents the position of an individual sensor 
as a function of time. The moving sensors are separated by 
a distance of A apart and move in the positive x direction 
at a constant speed of v represented by the slope tan6> of 
the lines in the figure. From ([25l) we know that the signal 
seen by each sensor is bandlimited to pt -\- vpx. We assume 
that the sensors sample in time at the temporal Nyquist rate 
of ^i±^. We further assume that the samples taken by the 
various sensors are all synchronized in time such that the 
collection of all samples lie on a two-dimensional lattice. In 
this case, we know from classical sampling theory | 6| that we 
can perfectly reconstruct any spatio-temporal field bandlimited 
to Q from its values at these sample locations provided that 
the repetitions of its spectra do not overlap in the spectrum of 
the samples. Furthermore, since the temporal sampling rate is 
above the Nyquist rate, it can be shown that only repetitions 
along one direction need be considered. These repetitions are 




(a) Sensor trajectories. 



-pt 



Pt 



> 



V 

(b) Repetitions in sampled spectrum. 
Fig. 6. Uniform configuration of moving sensors and resultant spectrum. 



illustrated in Figure |6(b)| for a rectangular set of the form 
^ = [—Px^Px] X [—Pt^Pt]- A detailed explanation of this no- 
alias condition can be found in |[T4ll for time-invariant fields. 

We now explicitly compute the no-alias conditions for two 
specific choices of Q. 

1 ) Spatio-temporal field bandlimited to rectangular region: 
Suppose ^} = [—px^ Px] X [—Pt^ Pt] is a rectangular region. Un- 
der the sampling configuration described above, the condition 
to ensure that there is no aliasing in the field reconstruction is 
that the repetitions shown in Figure [6(b)| do not overlap. This 
means that A should satisfy either 

27r ^ 27r , ^ 

— > 2px or —tdin0>2pt. 

Since tan^ represents the velocity of the sensors, the above 
condition is equivalent to the following requirement on the 
spatial separation between adjacent moving sensors: 



A < TT max <{ — , — /> . 

.pt Px 



(26) 



Thus Trmax is the maximum admissible spatial 

separation between adjacent sensors. 

2) Wave field: We now consider a time- varying field with 
a non-rectangular frequency spectrum. Suppose we are inter- 
ested in reconstructing the spatio-temporal wave field along a 
line. Assume that the field is produced by bandlimited sources 
located far from the region of interest. In this setting, we can 
use the far-field approximation to study the spectrum of the 
bandlimited field. It was shown in |24| that the spectrum of 
such a field is approximately supported on the region shown in 
Figure [7(a)| Here pt is the bandwidth of the source signals and 
Px = ^ where c denotes the speed of propagation of the wave. 
Now suppose that we sample the field using moving sensors 
with trajectories shown in Figure |6(a)| under the sampling 
configuration described before. Then, as before the condition 



10 




(a) Wave field spectrum (b) Sampled wave field spectrum 



Fig. 7. Far- field spectrum of a bandlimited source and i ts sampled version 
from samples taken by the mobile sensors of Figure |6(a)| 

required to ensure that there is no aUasing is that the spectral 
repetitions in Figure |7(b)| do not overlap. It follows from the 
figure that for sensor velocities v < c a. sufficient condition 

on the spacing A to ensure that the spectral repetitions do not 
27r px 2ii . t 1 

overlap is that — — < —2px^ r- tan 6', or equivalently 

A pt /\ 

TT V 

A <-(! + -). (27) 

Px C 

A similar analysis can also be performed for sampling a 
wave field over two-dimensional space using moving arrays 
of sensors. The sampling of such a field using an array of 
sensors is described in [251. These ideas can be extended to 
the case of mobile sampling by following an approach like the 
one we described in this section. 

B. Comparison with static sampling 

In Section HFdI we noted the advantages of mobile sampling 
over static sampling obtained by using an anti-aliasing filter to 
limit the contribution of out-of-band noise while sampling and 
reconstructing time-invariant fields. For sampling time- varying 
fields, however, such an advantage is not as significant. Con- 
sider filtering and sampling a one-dimensional time-varying 
field using sensors moving according to the configuration 
depicted in Figure [6(a)] In this case, we are essentially filtering 
along the lines in the t-x plane shown in Figure |6(a)| We note 
that it is possible to filter over time even with static sensors. 
This would amount to filtering along lines parallel to the t- 
axis in Figure |6(a)[ Thus the only difference between filtering 
in the mobile and static sensing cases is in the direction of 
filtering in the t-x plane. Hence the relative advantages of the 
two schemes would depend on the spectral characteristics of 
the additive noise. However, mobile sensing offers a different 
sort of advantage over static sensing: In some situations, we 
can get a reduction in the spatial density of sensor deployment 
required while using mobile sensors, as we show below. 

Consider a time- varying field bandlimited to a rectangular 
region [—px^px] x [—pt^pt] as in Section ITlI-All We know 
from classical sampling [6| that for sampling with static sen- 
sors the maximum spacing allowed between adjacent sensors 
is ^ . Comparing with the spacing requirement in the mobile 
setting given in ([26l) . it follows that when the mobile sensors 
are moving at a speed v > the inter-sensor spacing can be 

increased by a factor of ^ . In other words, this means that for 

^ pt 

a given length of the spatial region of interest, we can reduce 




5.33 5.3305 5.331 5.3315 5.332 5.3325 
eastings (m) 



Fig. 8. Spatial temperature field on EPFL campus. 



nSv/h 




Fig. 9. Bandlimited approximation of the radiation field near Fukushima. 

the number of sensors required by a factor ^ . The advantage 
is more significant when v ^ i.e., for spatio-temporal 
fields that vary slowly in time and at fast rate over space. 
This matches with the intuition that slowly varying fields are 
easier to track using mobile sensors. However, as the speed 
V is increased, the required temporal rate of sampling given 
by the Nyquist rate, also increases. Hence, in short, 

by using moving sensors we can reduce the spatial density of 
sensors at the cost of increasing their temporal sampling rates. 

Now consider the scenario of sampling a wave field along a 
line located in the far field of bandlimited sources as in Section 
IIII-A2I For sampling with static sensors the maximum spacing 
allowed between adjacent sensors is again Hence it follows 
from (|27l) that the inter- sensor spacing can be increased by a 
factor of (1 + ^) when we employ sensors moving at speed 
V. This observation suggests that for wave field reconstruction 
we get a significant improvement in the sensor spacing using 
mobile sensing only when the sensors can move at a speed of 
the order of the speed of wave propagation in the medium. 

IV. Simulations: Merits of mobile sensing 

We simulated the static and mobile sampling schemes for 
measuring the surface temperature field on a portion of the 



TABLE II 

Percentage root-mean-square errors with various schemes. 



Data type 


Static 
sensing 


Mobile 
sensing 


Mobile sensing with 
oversampling 


Temperature 


0.53% 


0.45% 


0.42% (no filter) 


Bandlimited radia- 
tion (SNR 20 dB) 


9.9% 


1.5% 


1.5% (with filter) 
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EPFL campus. For the true temperature field we used the 
readings obtained from ll26l as illustrated in Figure [H For 
static sensing, we considered sensors on a rectangular grid. 
For mobile sensing we assumed that the sensors move parallel 
to the X-axis and apply an anti-aliasing filter prior to sampling, 
like in Section III-B2I They take samples at the same points 
on the rectangular grid as in the static case. As seen in Figure 
[8] the field has sharp variations in space and hence is not 
bandlimited. Thus we expect some aliasing in the reconstruc- 
tion obtained via sine interpolation from samples of such 
a field. In Table Ull we list the percentage root-mean- square 
errors in the reconstructed fields, defined as x 100. 

The last column represents the performance obtained with 
mobile sensing assuming that the sensors measure the field 
at all points on their paths without any filtering. As the 
values in the table indicate, mobile sampling outperforms static 
sampling, and oversampling along the trajectories improves 
the performance further. The filtering operation in the mobile 
sampling scheme reduces the amount of aliasing in the samples 
leading to a reduction in the reconstruction error. We note 
that the temperature field is not truly bandlimited and hence 
the performance gains are more modest than what could be 
expected if the field were truly bandlimited. 

We also simulated the same schemes for sampling and 
reconstructing a truly bandlimited field in noise. For the true 
field shown in Figure [9] we used a bandlimited approximation 
to the spatial radiation field around the site of the Fukushima 
nuclear accident on 11 March 2011. The radiation levels in 
this field were measured in units of micro- Sieverts per hour 
(fiSw/h) at various positions during the months of July - 
September, 2011, and are available online at ll27ll . We consid- 
ered the sampling of a noisy version of this field with the noise 
spectrum as described in the statement of Proposition 12. II with 
the ratio of the sides a = 40 and estimated the percentage root- 



mean- square errors 



VE[||/-/||i] 



X 100 in the reconstruction. 



The values of the errors shown in Table HIl suggest that the 
reduction in the error obtained with mobile sensing is more 
significant than that was seen for the temperature field. We 
also see that the ratio of the errors under the static and mobile 
reconstruction schemes is approximately ai as expected by the 
result of Proposition 12.11 In this example we allow filtering 
in the oversampling scheme since the field of interest is 
bandlimited. From the last column of Table HIl we see that there 
is no improvement in accuracy with oversampling. This is 
expected since there is no advantage in increasing the sampling 
rate beyond the Nyquist rate. 

V. Conclusion and future work 

In this paper we have studied strategies for sampling and re- 
constructing a bandlimited spatial field using moving sensors, 
including both time-invariant and time- varying fields. We high- 
lighted and quantified the advantages of mobile sensing over 
classical static sensing, both in theory and through simulations. 
Our results for time-invariant fields clearly demonstrate the 
following advantages when using a time domain anti-aliasing 
filter together with a mobile sensor: 



(i) For non-bandlimited fields: Anti-aliasing filtering with 
mobile sampling suppresses aliasing in the direction 
of motion. For one-dimensional fields, higher sampling 
rates yield lower distortion in the reconstruction. 

(ii) For bandlimited fields in noise: Anti-aliasing filtering 
with mobile sampling suppresses noise in the direction of 
motion. This prevents aliasing in the direction of motion. 
Sampling at the Nyquist rate is sufficient. 

In the latter case we quantified the SNR improvement of the 
mobile sampling scheme over the static sampling scheme. 
For time- varying fields, we demonstrated the improvement in 
sampling density that can be obtained by using mobile sensors. 

In our analysis of mobile sensing of time-invariant fields 
in M? we considered only trajectories composed of a set 
of equispaced parallel lines. The analysis of SNR in the 
reconstruction can be generalized to higher dimensional spaces 
and to more general configurations of straight line trajectories 
like the ones studied in |[T4ll . In practice one may be forced 
to use non-linear sensor trajectories. In such cases, a possible 
approach would be to approximate the paths by straight line 
segments and then use the appropriate anti-aliasing filters to 
reduce spatial anti-aliasing as we did for the non-uniform 

However, quantifying 



speed sensing example of Section III-FI 
the noise suppression in such cases would be more complex. 
For general trajectories, it would be interesting to study the 
tradeoff between the SNR improvement and the path density 
metric of the trajectories introduced in |[T4ll . Similar extensions 
are also relevant for studying time- varying fields. 

For making mobile sampling schemes practical, one would 
also need to consider the effects of mobility on the field 
and on the sensing process. It is possible that the physical 
process of moving the sensor through the field may affect the 
characteristics of the field or introduce noise and irregularities 
in the sensing process. These effects must also be taken into 
account to completely characterize the advantages of mobile 
sensing over static sensing in practice. 
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Appendix 

Bandwidth of time-warped signal soit) = f{x{t)) 
Consider a piecewise-affine function of the form 

K 

Mt) = ^{uk^Vkt)X{t e [tk^tk^i)}. (28) 

k=l 

where tk < tk+i, and ^ < Vk < v where v denotes 
the maximum speed of the sensor. Let A/e := tk-\-i — tk, 
A := min/e Ak and tk := We know that the Fourier 

transform of si{t) := f{xi{t)) is given by 

K 



k=l 



1 2f±l^ f ^ 
Vk \Vk 



sine 



27r 



(29) 
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where F{.) denotes the Fourier transform of the field /(.), and 
the operation *^ denotes convolution with respect to ^. From 
the structure of the Fourier transform we can argue as in 1281 
that the effective bandwidth of 5i is given by 

1 _ 1 

max[v/eP + 7—] <vp^—. (30) 

k ZA/g Za 

We now study how the Fourier transform of the observed 
signal gets modified by a slight deviation from the piecewise 
affine trajectory. Suppose the trajectory is given by 

K 

= ^{uk ^ Vkt ^ eXk{t))X{t e (31) 

k=l 

Let px denote the maximum low-pass bandwidth of all Xk. 
Assume that / is twice differentiable. Then we have by 
Taylor's approximation that for t G [t/c,t/c+i), 

so{t) = f{x{t)) 

= f{uk^vkt^exk{t)), for t e[tk,tk+i) 
= f(uk + Vkt) + f\uk + Vkt)exk(t) + O(e^) 

Using Sq{£) to denote the spectrum of 5o(t) it follows that 
for small e we have 

e^A,sinc(^)|+0(e2).(32) 

Thus, we can argue that for small e the difference between 
the Fourier transforms of so(t) and si{t) is given by a term 
proportional to e over frequencies in the range 

1^1 < max[v/ep/ + + 7^] < ^P/ + P5 + X 

k ZA/e Za 

and only by terms of order O(e^) for other frequencies. By 
considering higher order terms in the Taylor series expansion, 
it follows that the difference between Fourier transforms of 
so(t) and si{t) is given by a term of order 0(e^+^) over 
frequencies outside of the range 

1^1 < vpf^mpx + ^. 
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